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Abstract 



Relaxation of electronic spins in metals is significantly enhanced whenever 
a Fermi surface crosses Brillouin zone boundaries, special symmetry points, or 
lines of accidental degeneracy. A realistic calculation shows that if aluminum 
had one valence electron, its spin relaxation would be slower by nearly two 
orders of magnitude. This not only solves a longstanding experimental puzzle, 
but also provides a way of tailoring spin dynamics of electrons in a conduction 
band. 
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Electronic spin is emerging as a building block of new digital devices. All-metal bipolar 
spin transistor |]J was already demonstrated and new devices are being built using the 
phenomenon of giant magnetoresistance ||. Recent advances in spin-coherent dynamics in 
semiconductors || may advance the development of a quantum computer Q]. All these new 
technological applications rely on relatively long relaxation time of conduction electron spin 
eigenstates. In this Letter we consider theoretically the problem of electronic spin relaxation 
in metals, and solve an outstanding puzzle, explaining why very similar metals (for example 
Al and Na) may have spin relaxation rates differing by two or three orders of magnitude. 

Experiments [||,^| show that spin states in metals live several orders of magnitude longer 
than momentum states. Furthermore, unlike momentum relaxation times r, spin relaxation 
times T\ [|7| vary significantly among metals. According to Elliott || and Yafet two 
factors cause electronic spin in metals to decay: (i) the spin-orbit interaction induced, for 
example, by (crystal) ions or impurities, and (ii) a momentum relaxation process such as 
impurity or phonon scattering. The product c 2 T 1; where c 2 is a measure of the strength of 
the spin-orbit interaction, is then well approximated by r and have similar magnitude for all 
simple metals. Considering only ion-induced spin-orbit interaction, crystalline and atomic 
c 2 should be similar. Indeed, by substituting the atomic values for c 2 Monod and Beuneu 
Tcfl found a "Gruneisen" behavior of c 2 T\ for all alkali and noble metals: as a function of 



reduced temperature T/T D , where T D is a Debye temperature, the values of c 2 Ti fall onto a 
single curve. On the other hand, metals Al, Pd, Mg, and Be have c 2 T x much smaller (up to 
three orders of magnitude for Mg and Be) than the "main group." For example, the atomic 
c 2 for aluminum differs by less than 10% from that for sodium [10|, yet the spin relaxation 



times at Debye temperatures are 0.1 ns for aluminum [|Tl|] (Tjj = 394 K |I2|| ) and 20 ns for 
sodium 0,|13|](Td = 150 K |12|]); the corresponding momentum relaxation times are in the 
ratio 1:7 |14|]. We resolve this puzzle by showing that the crystalline c 2 is about thirty times 
greater in Al than in Na due to a rather subtle "band renormalization" effect. 

In this paper we answer the question: "Why spin in some metals decays unexpectedly 
fast?" by introducing the concept of band renormalized spin-orbit interaction strength c 2 . 



Common to the "strange" metals (Al, Pd, Mg, and Be) with unusually fast spin relaxation is 
that their Fermi surfaces contain regions where c 2 is significantly enhanced. Such regions are 
found near Brillouin zone (BZ) boundaries, special symmetry points, or lines of accidental 
degeneracy. Although these spin "hot spots" comprise, generally, only a small part of the 
Fermi surface, they almost entirely determine the effective value of c 2 . We predict that the 
"strange" behavior is, in fact, common to all polyvalent metals, where spin hot spots are a 
consequence of the Fermi surface topology. This prediction is particularly significant since 
no polyvalent metals other than Al, Pd, Mg, and Be have been measured for 7\ so far. 
Furthermore, our calculations show that T\ can be appreciably altered by modifying the 
band structure (by alloying, doping, reducing the dimensionality, etc.). Our prediction is 
based on a realistic pseudopotential calculation for aluminum and analytical estimates for 
c 2 . 

If the periodic potential due to ions in a crystal lattice contains spin-orbit coupling 
(a term proportional to the scalar product of the orbital and spin momentum operators, 
L • S), the electronic Bloch states are a mixture of spin up | j) and down |j) species |5): 
*L»(*) = Kn(r) |T)+& k n(r) ||)]exp(ik-r) and tfjjr) = [a*_ kn (r) \l)-b*_ kn (v) |T>] exp(*k-r). 
The lattice momentum k is confined to the first BZ, n is a band index, and a^ n (r) and b kn (r) 
are complex periodic functions with the period of the lattice: if G denote the reciprocal 
lattice vectors, then ak n (r) = J2g a kn(G) exp(zG • r) and similarly for b kn (r). Both states 
have the same energy E^, as follows from time and space inversion symmetry ||; the 
numbering of bands is therefore the same as without the spin notation. The degenerate 
states and \I/ kn are chosen to represent electrons with spins polarized along z direction 
0: (^knl^l^kn) = — (^Lil^l^kn) < and the off-diagonal matrix elements are zero. 
This condition implies that a kn (r) have values close to one, while &kn( r ) are much smaller, 
decreasing with the decrease of the strength of the spin-orbit interaction (with the exception 
of the points where the spin-orbit interaction lifts a degeneracy). 

Elliott H noticed that ordinary (spin conserving) impurity or phonon scattering can 
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induce transitions between W kn and ^j^/, leading to the flip of a spin polarization and 
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thus spin relaxation. If (b 2 ) is the Fermi surface average of |&kn| 2 = J2g |^kn(G)| 2 , spin 
relaxation rate 1/Ti « 4(6 2 )l/r (that is, c 2 ~ A(b 2 )). Such a formula could be obtained 
by assuming that (1) 1/r ^> 1/Ti (2) akn(r) ~ 1, (3) scattering form factor has a 
constant amplitude, and (4) the interference between &k«(G) with different G is neglected. 
Assumptions (1) and (2) are usually satisfied. Assumption (3) implies a scattering by a 
delta-function-like impurity potential. This simplification is consistent with our goal to 
establish the effect of band-structure (through b^ n ) on Ti rather than to study particular 
scattering processes. The neglect of interference can be justified if the scattering form factors 
have rapidly varying phase. Although such a form factor [that also satisfies (3)] is hardly 
found, realistic form factors do oscillate on the momentum scale where the transitions occur. 
Assumption (4) is then partly justified, since the phase attached to &kn(G) will be different 
for different G. 

The spin-mixing parameters |&kn| 2 can have a broad range of values, depending on the po- 
sition of k in the BZ. Consider a band structure computed without the spin-orbit interaction. 
If the closest band to n is separated from n by A, the spin-orbit interaction mixes the spins 
from the two bands (direct interband transitions): |&kn| 2 ~ (1 — A/y^ 2 + 4V| )/2, where 
Vso is some effective spin-orbit interaction. Three cases occur. (A) For a general point, the 
band separation is of order Ep, the Fermi energy, so that A ^> Vso an d |&k«| 2 ~ (Vso/Ep) 2 . 
(B) If the state is close to a BZ plane that cuts G by half, the band separation is ~ 2Vg (Vq 
is the Gth Fourier coefficient of the non-spin part of the lattice potential). Since typically 
Vb ^> Vso, \bkn\ 2 ~ (Vso/2Vg) 2 ; this can be a few orders larger than in (A). Finally, (C) 
the spin-orbit interaction can lift the degeneracy of two or more bands. The mixing of spins 
is complete and |6kn| 2 ~ |«kn| 2 ~ 0.5. The states with property (C) are not covered by the 
Elliott- Yafet theory 0. This is not a concern, however, since as we show below, such states 
are statistically irrelevant for aluminum Ti. 

To illustrate how the band structure affects (b 2 ) (and Ti), we perform a pseudopotential 
calculation for aluminum, where all three cases (A) to (C) occur. Our pseudopotential has 
a non-spin (scalar) part and a spin-orbit term. The scalar part [|l6j has the nice feature [in 
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view of estimate (B)] of being fitted to the experimental values [17]] of V\ = Vni = 0.00895 
a. u. and V 2 = V 20 o = 0.0281 a. u. (1 a. u. = 2 Ry). In addition to accurately reproducing the 
band structure of aluminum this form factor gives reasonable results ||16[ for the resistivity 
of liquid aluminum, making it useful in scattering problems. The spin-orbit part of our 
pseudopotential is AL • SVi, where Vi is the operator projecting on the orbital momentum 
state /. The parameter A = 2.7 x 10~ 3 a. u. inside the ion core of twice the Bohr radius, 
r c = 2tb- Outside the core A vanishes. This estimate for A is based on the first-principles 



spin-orbit pseudopotential of Ref. ||18|| . Since the spin-orbit interaction acts only inside 
the core, the effective Vso is reduced to about (r^/V a )X (~ A/10 for aluminum, where the 
primitive cell volume V a ~ lllr^). 

Having the band structure, the challenge is to evaluate the average of |&kn| 2 , a quantity 
that varies over several orders of magnitudes. We do the averaging by the tetrahedron 



method [|19j with a carefully designed grid that envelops the two sheets of the aluminum 
Fermi surface (bands 2 and 3). The grid is denser in the regions where A is smaller (around 
BZ boundaries and accidental degeneracy points). This adaptation of the grid is necessary 
to ensure that |&kn| 2 can be linearly interpolated inside the grid cells (tetrahedrons), as 
assumed in the tetrahedron method. The number of grid points (in the irreducible wedge 
of the first BZ) where the band-structure equations must be solved is about ten thousand. 

Figure [l] plots the calculated distribution p of the values attainable by |&kn| 2 over the 
Fermi surface (the visual perspective of the distribution is in Fig. |2|). The span is enormous- 
almost seven decades! The majority of states have |&kn| 2 below 10~ 5 . These are the generic 
points from the estimate (A). Once the Fermi surface approaches the BZ planes (violet in 
Fig. |2|), the values jump to 10~ 5 — 10~ 4 . However, the largest |&kn| 2 are found near the 
accidental degeneracy points R (red spots in Fig. 0); these points do not lie on symmetry 
lines, yet they have degenerate bands []20| if no spin-orbit interaction is present. The spin- 



orbit interaction lifts this degeneracy and completely mixes the spin states as in the estimate 
(C). The unusually long tail of p ensures the large value of the average (b 2 ) ~ 2.0 x 10~ 5 , 
about ten times larger than the value where p is maximal. To ensure that the above picture is 



valid for a large temperature range, we checked that spin hot spots survive energy excitations 
of at least ~ kTp> (« 35 meV) above and below the Fermi surface. This also invalidates the 
objection that small relativistic corrections (of a few meV) could affect the existence of the 
spin hot spots in aluminum. 

The band structure role in (b 2 ) becomes even more evident from the plot (Fig. [I]) of 
p for a hypothetical case of monovalent aluminum (the lattice and the form factors are 
unchanged). The monovalent aluminum has a simple Fermi surface, without deformations 
of type (B) or (C). The distribution of b 2 is appropriately narrower, and the average value 
(b 2 ) ~ 3.4 x 10 -7 is about fifty times smaller than for trivalent aluminum. Adjusting for 
the density of states (monovalent aluminum would have 1/r reduced ~ 1/3 1 / 3 times) the 
spin relaxation would be about seventy times slower. These values may somewhat vary for 
different scattering processes. 

How different spin "hot spots" contribute to the renormalization of {b 2 ) becomes clear 
from Fig. ^, which plots what we call the average density of bands (ADOB). ADOB is 
the Fermi surface average of the number of (vertical) bands in the interval (A, A + dA): 
ADOB (A) = (l/g F ) Ek^m^kn - E F )5(\E km - E F \ - A), were band energies E kn are 
computed without the spin-orbit interaction and g F is the density of states (per spin) at 
E F . The average (b 2 ) can then be well approximated by taking the integral (from zero to 



infinity) of ADOB weighted by the mixing factor (1 — A/yA 2 + 4V| )/2. Figure § uses 
Vso — 11 X 10~ 4 a. u., the value that gives the right answer for (b 2 ). At small A ADOB 
is linear. This is expected for a region around R where band gaps increase linearly with 
increasing distance from R |2(J. The linear increase (see also Tab. |) continues up to the 
point where the Fermi surface crosses the first (A « 2Vj) and the second (A « 2T4) closest 
BZ boundary plane. At these points ADOB has power-law singularities (Tab. |). Larger 
band separations, where ADOB develops either logarithmic or step-like singularities pl| , are 



irrelevant. Indeed, the cumulative average (b 2 )(A) saturates after the second peak so that 
(6 2 ) is almost entirely determined by the regions close to accidental degeneracies and BZ 
boundaries. As Fig. ^ shows, these regions contribute about equally to (b 2 ). The behavior of 



(6 2 )(A) also shows that states in the immediate neighborhood of R (the red spots in Fig.0) 
with A < 2Vso ~ 2.2 x 10~ 4 are statistically irrelevant. It is rather a broader neighborhood 
of R [the states with 2Vso ~ A < 2V\ and estimate (B)] that are contributing most. These 
findings are confirmed by the analytical calculation reported in Tab. |I[ 

Table [I] summarizes our estimates of AD OB and (b 2 ) for three different cases. The 
estimates were obtained analytically by the orthogonalized-plane-wave methods developed 
in Ref. ||20|| . If the Fermi surface (in an extended zone) lies entirely within the first BZ (region 
I, alkali and noble metals), the enhancement of {b 2 ) is about tenfold since Emin ~ 0.1. If 
the Fermi surface crosses a zone boundary (region II, polyvalent metals) the enhancement is 
~ 1/Vq, typically a hundred. Relative to I, however, the enhancement is in the units often, 
in agreement with the numerical calculation. Qualitatively, a fraction of ~ Vg states on the 
Fermi surface comes close to zone boundaries (A ps 2V"g) so that their |&kn| 2 ~ (Vso/Vg) 2 - 
The contribution of these points to (b 2 ) is therefore ~ Vg x (Vso/Vq) 2 ] the enhancement is 
~ 1 /Vq, consistent with the result in Tab. |. Curiously, the region II applies also to the noble 
metals, whose Fermi surfaces touch zone boundaries ]I2" |. This effect is, however, masked 
by the unusually high Vg for these metals [^2J (Vg ~ Emin, otherwise the Fermi surface 
would not touch the zone planes). The enhancement of (b 2 ) can be also significant if the 
Fermi surface crosses a line of accidental degeneracy (region III) . Indeed, by substituting the 
values for aluminum, the enhancement is about tenfold, similar to case (II), again agreeing 
with the numerical result. Finally, not shown in the Tab. H is the case when Ep coincides 
with a degenerate level at a special symmetry point. Such a situation occurs, for example, 
in Pd and Pt, whose Fermi surfaces go through the fee L point. If the spin-orbit interaction 
lifts this degeneracy, the renormalization of (b 2 ) can be significant (we find the enhancement 



Vg/Vso for the fee W point Wl 



7 



TABLES 



FS 


multiplicative 


ADOB 




area 


factor 






I. 


N G /AGk F 


0(A - Emin) 


1/Emin 


II. 


N G /AGk F 


A/^A 2 - AV 2 


(7r/4)(i/y G ) 


III. 


N R /8k F 


(L 2 /Vi) 2 A 


{V2/Vif\n(Vi/V S o) 



TABLE I. Estimated contributions of different Fermi surface regions to ADOB and (b 2 ). (I) The 
Fermi surface is entirely within the first BZ, crosses (II) a zone boundary, and (III) an accidental 
degeneracy line at R. Momentum and energy is measured in the units of G/2 and (G/2) 2 /2, 
respectively. The region III assumes an fee lattice with period a, energy and momentum in the 
units of 2n/a and (27r/a) 2 /2, respectively, and L 2 > V\ > V 2 . Both ADOB and (b 2 ) come with the 
corresponding multiplicative factor. Notation: k F = \JE F is the Fermi vector, Emin = G{G—2k F ) 
is the band gap at the point of closest approach to the BZ plane given by G; Nq, and Nr are the 
numbers of the corresponding BZ planes and accidental degeneracy points (A^ = 24); O is the 
step function. 
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Our final note concerns the hexagonal Mg and Be, where the deviation of T\ from the 



"main group" is most striking |I0| . We argue that this is also a manifestation of the band 
renormalization of c 2 . Without the spin-orbit interaction, all the states at the hexagonal 
faces of the first BZ of a simple hexagonal structure are degenerate ||12| . The spin-orbit 
interaction lifts this degeneracy || (except at some symmetry points and lines), presumably 
by the amount VgVso, the largest second-order term containing the spin-orbit interaction 
(any first order term vanishes since the structure factor associated with the hexagonal faces 
is zero [|12|]). The contribution to (b 2 ) of the points where the Fermi surface intersects the 
hexagonal faces is ~ VgVso (in the units of Tab. |): the characteristic value |6kn| 2 ~ 1, 
times the area of the affected part of the Fermi surface, VgVso- The enhancement measured 
in terms of Vg Q is then Vg/Vso', this can be as large as a thousand for light elements like 
Mg and Be. 

We acknowledge discussions with P. B. Allen and M. Johnson. This work was supported 
by the U.S. ONR. 

Note Added.-Aiter submission of our work R. H. Silsbee brought to our attention an 
earlier paper (R. H. Silsbee and F. Beuneu, Phys. Rev. B 27, 2682, 1983) which suggests 
the importance of the accidental degeneracy points for the spin relaxation in aluminum. 
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FIGURES 

FIG. 1. Calculated distribution p (in arbitrary units) of the spin- mixing parameters |&kn| 2 for 
aluminum. The corresponding average (b 2 ) ~ 2.0 x 10~ 5 is indicated by solid arrow. The linear 
tail of the distribution is shown in the inset. The dashed line shows what the distribution would 
be if aluminum were monovalent ((b 2 ) pa 3.4 x 10~ 7 , dashed arrow). 

FIG. 2. Stereogram of the Fermi momentum directions in aluminum. The fragment shows the 
spin "hot spots": the points k (in an extended-zone scheme) with |6kn| 2 > 10~ 5 . Colors violet, 
blue, green, yellow, and red indicate a successive increase of b 2 : violet points have b 2 from 1CT 5 to 
1CT 4 , blue 10~ 4 to 10~ 3 , etc., up to 10" 1 — 1 for red. To improve their visibility, the weight of the 
colors (except for violet) is enhanced. 

FIG. 3. Calculated average density of bands (ADOB) and cumulative average spin-mixing 
parameter (b 2 )(A) for aluminum. 
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